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O . ABSTRACT 



Context. Cometary globule (CG) 1 and CG 2 are "classic" cometary globules in the Gum Nebula. They have compact heads and long 
dusty tails that point away from the centre of the Gum Nebula. 

Aims. We study the structure of CG 1 and CG 2 and the star formation in them to find clues to the CG formation mechanism. The 

■ two possible CG formation mechanisms, radiation-driven implosion (RDI) and a supernova blast wave, produce a characteristic mass 
' distribution where the major part of the mass is situated in either the head (RDI) or the tail (supernova blast). 

Methods. CG 1 and CG 2 were imaged in the near infrared (NIR) JsHKs bands. NIR photometry was used to locate NIR excess 
objects and to create visual extinction maps of the CGs. The Ay maps allow us to analyse the large-scale structure of CG 1 and 

■ CG 2. Archival images from the WISE and Spitzer satellites and HIRES-processed IRAS images were used to study the globule's 
small-scale structure. Fits were made to the spectral energy distribution plots of the NIR-excess stars to estimate their age and mass. 
Results. In addition to the previously known CG 1 IRS 1 we discovered three new NIR-excess objects in IR imaging, two in CG 1 

^ • and one in CG 2. CG 2 IRS 1 is the first detection of star formation in CG 2. The objects are young low-mass stars. CG 1 IRS 1 is 

■ probably a class I protostar in the head of CG 1. CG 1 IRS 1 drives a bipolar outflow, which is very weak in CO, but the cavity walls 
^ I are seen in reflected light in our NIR and in the Spitzer 3.6 and 4.5 yum images. Strong emission from excited polycyclic aromatic 

hydrocarbon particles and very small grains were detected in the CG 1 tail. The total mass of CG 1 in the observed area is 41.9 Mq 
of which 16.8 Mg lies in the head. For CG 2 these values are 31.0 Mq total and 19.1 Mgin the head. The observed mass distribution 
does not offer a firm conclusion for the formation mechanism of the two CGs: CG 1 is in too evolved a state, and in CG 2 part of the 
^ ' globule tail was outside the observed area. 

CO ■ Conclusions. Even though the masses of the two CGs are similar, star formation has been more efficient in CG 1. By now, altogether 

■ six young low-mass stars have been detected in CG 1 and only one in CG 2. A possible new outflow was discovered to be emanating 
CS| ■ fromCGlIRSl. 

Key words. Stars: formation - Stars: pre-main-sequence - ISM: individual (CG 1, CG 2) - ISM: dust, extinction 

1. Introduction ture. Determining the triggering mechanisms leads to a better 

understanding of star formation processes and the characteris- 

Star formation takes place in dense molecular clouds when grav- ^j^^ newly formed stars 

ity starts dominating over the gas pressure. In triggered star for- ^ , , , , • , , ,- , , , 

, . mation this first "kick" to create a density fluctuation that will ^^^^^^^^^^^ globules (CGs) are a special class of globules. 

H . lead to the cloud collapse can be provided by mechanisms such ^Gs have dusty compact heads and elongated, faindy lumi- 
• as stellar winds or supernova explosions (e.g. Elmegreen 1998). ^f' ■ ^hey were fii-st discovered m the Gum Nebula 

The coflapse continues until the conditions in the coflapsing (Hawarden & Brand 1976; Sandqvist 1976 . CGs are typically 

cloud are suitable for star formation and the gas pressure starts ^""""f^ associations in HII regions (e.g^Rosette, White 

to resist further collapse. Small isolated molecular clouds and 1997 ■ Some of the classical CGs in the Gum Nebula are 

globules, such as the ones described by Bok & ReiUy (1947), ^^^es of isolated low-mass star formation (e.g. Reipurth 1983). 
are ideal locations to study the star formation process in detail. Two formation mechanisms for CGs have been suggested. 

The advance of large-format NIR arrays and especiafly of space- Brand et al. (1983) argued that CGs form when a supernova 

borne telescopes like the Spitzer Space Telescope has made it (SN) explosion shocks an originally spherical cloud. The shock 

possible to detect deeply embedded low-mass young stellar ob- compresses the cloud to form the head, and the blast wave 

jects in seemingly starless dark clouds (see e.g Young et al. drives the material mechanically away from the SN to form the 

2004). These embedded objects are not visible at shorter wave- tail. Reipurth (1983) suggested radiation driven implosion (RDI) 

lengths because of the high extinction. If the formation of these where the UV radiation from massive O stars photoionizes the 

stars was triggered by some outside event, the triggering mecha- cloud and shock fronts compress the cloud. The less dense cloud 

nism may have left a telltale sign into the large-scale cloud struc- medium is separated from the cloud by radiation and ionization 

shocks, and the tail is formed out of the eroded cloud medium. 

* Based partly on observations done at the European Southern Numerical simulations have been done for both mechanisms and 

Observatory, La Silla, Chile (ESO programme 078.C-0490). they show that both RDI and SN can create tails behind the 

** Appendices A and B are only available in electronic form at globule core. According to the RDI model, a major part of the 

http://www.aanda.org CG mass is in the "head" (e.g. Kessel-Deynet & Burkert 2003; 
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Lefloch & Lazareff 1994; Bertoldi & McKee 1990), whereas in 
the SN model a large part of the mass can he in the tail (e.g. 
Heathcote & Brand 1983; Bedogni & Woodward 1990). This 
offers the possibility of determining the mechanism that triggers 
the formation of CGs by studying the mass distribution in them. 

We used NIR imaging to study the structure of and star for- 
mation in two archetype cometary globules, CG 1 and CG 2, 
in the Gum Nebula. These globules are a part of the system of 
cometary globules in the Gum Nebula (Reynolds 1976a). This 
nebula has been described as a supernova remnant and/or an 
HII region. A SN explosion ~1 Myr ago has been suggested as 
the possible origin of the Gum Nebula (e.g. Reynolds 1976b; 
Brand et al. 1983). As an HII region, the Gum Nebula is mainly 
ionized by two O stars, the multiple system Velorum and f 
Puppis. The tails of the CGs point away from the Gum Nebula 
centre, suggesting a common triggering mechanism (Hawarden 
& Brand 1976; Zealey et al. 1983). CG 1 is one of the largest 
of the Gum Nebula CGs. Its tail is 25' long and the characteris- 
tic size of the globule head is 2' (Brand et al. 1983). CG 2 is a 
smaller CG situated -52' northwest (NW) of CG 1. Even though 
the tail length is listed as 26' by Zealey et al. (1983) and 18' by 
Sridharan (1992), the effective tail length is and the rest is 
only marginally detected in the optical. The distance estimates to 
CG 1 range from 300 to 500 pc, and this paper adopts the value 
300 pc of Franco (1990) used in Haikala et al. (2010) (Paper I). 

The observations of the CG 1 head have been discussed in 
Paper I where the detection of a young stellar object embed- 
ded in the dense globule head was reported. This paper expands 
the observations to cover the CG 1 tail and the nearby CG 2. 
Observations and data reduction are described in Sect. 2, and the 
results are presented in Sect. 3. Discussion is in Sect. 4 and a 
sunnmary in Sect. 5. 

2. Observations and data reduction 

2.1. NIR Imaging 

CG 1 was imaged in J, H, and Ks with the Simultaneous 
InfraRed Imager for Unbiased Survey (SIRIUS, Nagayama et al. 
2003) on the InfraRed Survey Facility (IRSF) 1 .4m telescope at 
SAAO in Jan. 2007. The SIRIUS field of view is 7:7x 7^7 and 
the pixel scale 0'.'453. The observations were carried out in the 
on-off mode to preserve the source surface brightness. The off- 
position was at 7''18'" -44°05', (J2000). Each on and off ob- 
servation consisted of ten individual jittered images. Depending 
on the stability of the weather 10 s or 15 s integration time was 
used for the individual images. The average seeing was r.'3-r.'4. 
Four overlapping fields were observed to cover the area where 
'^CO was observed by Harju et al. (1990). Going from east to 
west, the fields will be referred to as 1(RSF)1, 12, 13, and 14. The 
observing sequence was off - II - 12 - off - 13 - 14 - off, and thus 
the total integration time (~1500 s) for each field is the same. 
Sky flats were observed in the evening. CG 2 was imaged in the 
same way as CG 1, but only two overlapping fields were ob- 
served. The total integration times of the fields are ~2200 s and 
-1200 s. 

CG 1 was imaged in /s, H, and Ks with SOFI (the Son 
of Isaac) near-infrared instrument on the New Technology 
Telescope (NTT) at the La Silla Observatory, Chile. The SOFI 
field of view is 4:92 X 4:92 and the pixel size is ^.'288. The 
imaging covered the area where Harju et al. (1990) detected 
'^CO emission. Three overlapping SOFI fields were observed. 
The fields will be referred to as N(TT)1, N2 and N3 (from east to 
west). The observations were carried out in Feb. 2007 except for 



the field N3 in Ks, which was observed in Jan. 2010. SOFI imag- 
ing was done in the standard jitter mode with jitter box widths of 
20"-30". Fields Nl and N2 were imaged in observation blocks 
of 50 individual images of one minute integration each. Three 
and two observing blocks in all colours were obtained for N 1 and 
N2, respectively. One observing block of 25 images in the three 
colours was obtained in field N3 and in the same off-field which 
was used in the IRSF observations. Standard stars from Persson 
et al. (1998) were observed before and after each observation 
block. CG 2 was imaged with SOFI in Js and Ks bands in Feb. 
2007 and in H band in Jan. 2010. The observations were done in 
the standard jitter mode with a jitter box width of 30". Images 
were obtained in observation blocks of 25 individual frames with 
one minute integration each. Js and H bands have one and Ks 
band has two observation blocks. The average seeing during the 
SOFI observations was 0'.'8. 

2.2. Data reduction 

The data for SIRIUS and SOFI was reduced as described in 
Paper I. For the sake of completeness, the data reduction steps 
for the single fields are included in Appendix A. Data from the 
separate frame catalogues were combined into one. Any dupli- 
cate stars from the overlap regions were averaged into a single 
entry. 

The final catalogues for SIRIUS fields contain 106 objects 
(11), 161 (12), 218 (13) and 172 (14). For the SKY field, the final 
catalogue has 313 objects. The CG 2 catalogues have 337 (II) 
and 186 objects (12). The final combined SIRIUS catalogue for 
CG 1 and CG 2 have a total of 576 and 466 stars, respectively. 

The number of objects in SOFI catalogues for CG 1 are 261 
in Nl, 296 in N2 and 142 in N3. The combined photometry cat- 
alogue of CG 1 has 674 stars, the single CG 2 frame has 201 
and the off-field 230 stars. The stellar limiting magnitudes cor- 
responding to a OT^l error in the SOFI catalogues for Nl and N2 
are 20?5 in H and Ks and IV^S in 7s. For N3, OFF, and CG 2 
the limiting magnitudes are ~0.5-l .0'" brighter. 

3. Results 

3.1. NIR Imaging 
3.1.1. CG 1 

A false-colour image of the SIRIUS data combined from on- 
off-mode observations is shown in the online Fig. B.l. The /, 
H, and Ks bands are coded in blue, green, and red, respec- 
tively. The brightest object in flie image is NX Pup at 7''19'°28.3' 
-44''35'H". A bright reflection nebula (RN) is seen below NX 
Puppis in the figure and is most probably illuminated by it. A 
thin, bright filament leads southwest (SW) from the RN and a 
fainter filament extends NW from the nebula. The reddish, semi- 
nebulous bright patch 70" west of NX Pup is due to reflected 
light from a newly born star (young stellar object; YSO) embed- 
ded in the dense globule head. Extended surface brightness is 
observed NW of the YSO. A bright arc reaches northeast (NE) 
from the YSO. See Paper I for discussion on the detailed struc- 
ture of the head. About 6' west of the YSO lies another region of 
extended surface brightness. The on-off observation mode pre- 
serves the true surface brightness and therefore the extended 
structure seen in the onUne Fig. B.l is real. 

A false-colour image of CG 1 combined from SOFI images 
is shown in Fig. 1. The objects CG 1 IRS 2 and 3 (see Sect. 
4.3) are indicated. NX Pup is too bright to be observed with 
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Fig. 1. Colour-coded SOFI image of CG 1. Tlie J, H, and A's bands are coded in blue, green, and red. The objects CGI IRS 2 and 3 (from east to 
west) are indicated with a box. 



SOFI and was therefore left outside the image. The edge of the 
RN below NX Pup is seen in the lower left-hand corner of the 
figure. The filamentary structure and the YSO are seen more 
clearly than in Fig. B.l. The jittering observing mode uses the 
observed ON frames to estimate the sky brightness. This smears 
out surface brightness structures larger than the jitter box width 
and produces filaments that trace gradients in the original back- 
ground. Point sources and galaxies are unaffected. Therefore the 
extended surface brightness NW of the YSO seen in the online 
Fig. B.l has largely disappeared, but the northern filament is 
seen more clearly because of strong local surface brightness vari- 
ation. The negative brightness areas right next to the filaments 
are not real but due to the jittering. The overall morphology of 
CG 1 in the SOFI data can be best seen in the Js and H bands 
(online Fig. B.2, upper panels). In Ks (online Fig. B.2, lower 
panel), only the leading edges, a short filament behind them, and 
the brightest tail filament are seen. Most of the extended sur- 
face brightness NW from the YSO seen in the SIRIUS image 
has been transformed into a faint filamentary structure. In the 
western part of the tail, another fainter, slightly arced filament is 



3.1.2. CG2 

A false-colour image of CG 2 obtained with SIRIUS in the on- 
off observing mode is shown in the online Fig. B.3. The J, H, 
and A's bands are coded in blue, green, and red, respectively (the 
individual SIRIUS J, H, and Ks images are shown in the on- 
Hne Fig. B.4). The object CG 2 IRS 1 (see Sect. 4.3) is indicated 
with a box in the image. In J (online Fig. B.4, top left panel), ex- 
tended surface brightness is observed in the globule head. There 
is a "nose" of surface brightness just east of a more extended 
region of enhanced surface brightness where visibly reddened 
stars are located. No details of the tail can be distinguished in 
the SIRIUS images. Only the NE part of the online Fig. B.3 is 
free of extended surface emission, and the sky to south and to 
east of the globule head is covered by faint surface brightness. 

The SOFI image of CG 2 is shown in Fig. 2 in false colour 
(the individual Js, H, and A's images are in the online Fig. B.5). 
The object CG 2 IRS 1 is marked. The "nose" appears slightly 
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Fig. 2. Colour-coded SOFI image of CG 2. The /s, H, and A^s bands are 
coded in blue, green, and red. The object CG 2 IRS 1 is indicated with 
a box. 



bluish. Behind the nose, a filamentary structure produced by the 
jittering observing mode extends to the NW. Reddened stars can 
be seen in this extended region. A faint filament runs along the 
SW edge of the globule. The CG 2 tail is not seen because it 
extends far outside the ~5' field of view. 



3.2. Photometry 

The SOFI limiting magnitude in each colour is ~2'" brighter than 
in the SIRIUS data. The SOFI data is discussed in detail in the 
following. Galaxies should have been removed from the cata- 
logue by the selection rules listed in the online Appendix A, but 
it is likely that some still remain. The catalogues containing the 
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SOFI photometry for CG 1, CG 2, and the OFF region are avail- 
able at the CDS. 

3.2.1. CG 1 

The J-H, H-Ks colour-colour diagrams of CG 1 and the OFF 
field are shown in the left and middle panels of Fig. 3. The 
Bessell & Brett (1988, BB88) main sequence converted into 
the 2MASS system using the conversion formulae by Carpenter 
(2001) is also plotted. The solid lines indicate the locations of 
unreddened main sequence. The dashed lines are reddening lines 
drawn using the BB88 reddening law; normal reddened MS stars 
lie between or close to these lines. The infrared excess stars 
should lie right of the rightmost reddening line. The discov- 
ered IR excess stars are named. The arrow indicates a reddening 
vector of 5?0. The errorbars indicate the average errors, which 
are about half of the maximum errors. Several stars are seen to 
the right of the reddening line in Fig. 3. In the ON field they 
are fainter than my ~ 20^2. Only one star brighter than this 
(J-H ^ 1.74, H-Ks= 1.24, m j =19?6) has infrared excess; 
this object (CG 1 IRS 2) is discussed further in Sect. 4.3.3. 

The OFF field shows no field-wide reddening effects, thus 
allowing its use as a reference field when studying the colour 
excesses. The objects to the right of the reddening line in the 
OFF field are faint (my > 197*3). The low signal-to-noise ratio 
of these objects (integration time of the sky field is lower than 
in the ON field) can cause Source Extractor (SB) to classify ex- 
tended objects as stars. The objects below the reddening Une are 
faint, and at a distance of 300 pc they would be substellar. It 
is possible that they are long-period variables or extragalactic. 
None of the objects showing infrared excess in the OFF field are 
detected in the Wide-field Infrared Survey (WISE) images (see 
Sect. 4.3 for more discussion on WISE). The objects with nega- 
tive (H-Ks) colours are possibly extragalactic. Both the ON and 
OFF field diagrams show a slight offset (H-Ks ~ 07*04) between 
the tabulated main sequence and the data. This offset is the same 
in both the ON and OFF regions, and thus it does not affect the 
reddening estimates. For comparison, 2MASS PSC data from a 
1000" box around the OFF position has a similar offset between 
the magnitudes and the tabulated main sequence. This suggests 
that the stellar colours differ from the tabulated colours at high 
gidactic latitudes, possibly because of metalUcity or because the 
transformation of the BBSS main sequence to the 2MASS pho- 
tometric system is not accurate. 

3.2.2. CG 2 

Except for the offset, the CG 2 J-H, H-Ks colour-colour dia- 
gram in the right panel of Fig. 3 shows good agreement with the 
main sequence and has also several reddened stars. The NIR ex- 
cess star at y - = 1 .73, - /:s = 1 .45 (CG 2 IRS 1) is located 
in the head of CG 2, just east of the globule core. The star is 
discussed in detail in Sect. 4.3.5. Stars below the reddening line 
are all fainter than 18"7 in nij, most of them fainter than 19?4. 

4. Discussion 

4.1. Visual extinction 

Visual extinction (Ay) maps were derived from the data using 
the NICER method introduced in Lombardi & Alves (200 1 ) . The 
method utilises colour excesses from the ON and the unreddened 
OFF field stars to determine the relative extinction between the 
two fields. 
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Fig. 4. Visual extinction map of CG 1. Upper panel: SIRIUS data. The 
contours are from 1 to 5 magnitudes in steps of 1'". Lower panel: SOFI 
data. The contours arc from 1 to 8 magnitudes in steps of 1". The grey 
scale is the same in both panels. 

The BBSS extinction law was used to estimate the visual 
extinction. The Mathis (1990) and Rieke & Lebofsky (1985) 
reddening laws produce ~9% and ~8% higher Ay values than 
BBSS, respectively. However, the relative variation in the ex- 
tinction within the final map will be the same irrespective of the 
reddening law used. The Ay values in the text have been derived 
with the BB88 law. 

In the SlRlUS data, the Gaussian used for smoothing has an 
FWHM of 72" and a pixel scale of 36", whereas the SOFI data 
was smoothed using a 45" FWHM Gaussian and the pixel scale 
is 22.5". The pixel values for any empty pixels were interpo- 
lated from the neighbouring pixels. Due to higher resolution and 
Umiting magnitudes, the extinction values in the SOFI map can 
reach higher values (and sharper peaks) than in the SIRIUS map. 

4.1.1. CG 1 

The extinction maps of CG 1 derived from the SIRIUS and SOFI 
data are shown in Fig. 4 upper and lower panels, respectively. 
The lowest contour value and the contour step in the figure is 
1?0. Typical values for the extinction outside the globule are less 
than l"" in the SIRIUS data. The highest extinction, S™3, is seen 
in the SOFI map. The large-scale structures in the SIRIUS and 
SOFI Ay maps agree well. Two large clumps are seen in both 
the SOFI and the SIRIUS maps and will be referred to as the 
"Head" and "Middle" clumps in the following. A further weaker 
clump called "Tail" is seen in the SIRIUS Ay map outside the 
area imaged with SOFI. The SOFI map has a small local maxi- 
mum at the western edge, but this cannot be seen in the SIRIUS 
map. This feature is most likely an artefact caused by the bright, 
saturated star at this location in the SOFI image. 

The Head maximum peaks west of the YSO and it agrees 
well with the location and shape of the maximum C'^O emission 
(see Paper 1) in the head of CG 1 . Towards north the eastern rim 
of the extinction in the head follows the sharp edge of the globule 
seen in Fig. 1. The Head maximum is thus located very close 
to the surface of the exposed side of the globule. The position 
of the Middle maximum agrees well with the extended surface 
emission seen in the SIRIUS image and the filamentary structure 
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Fig. 3. SOFI J-H, H-Ks colour-colour diagram of CG 1 (left panel), the OFF region (middle) and CG 2 (right). The dashed lines show the 
effect of the Bessell & Brett (1988) (BB88) reddening law. The unreddened main sequence given in BB88 is shown with continuous line. A 5™ 
reddening vector has also been included in the diagram. The average formal error of a data point is indicated in the lower right of the left panel. 
The discovered IR excess stars have been marked. The cross shows CG 1 IRS 3, which is not included in the final photometry catalogue. 



in the SOFI image. There is indication in the SOFI data that the 
Middle clump is fragmented into two. 

Unlike to what is seen in SOFT data, the maximum Ay in 
CG 1 in the SIRIUS data is not seen in the Head but in the 
Middle maximum. This happens because of the small size of the 
Ay maximum in the Head is diluted in the large SIRIUS pixels. 
In addition, because SIRIUS has lower limiting magnitudes than 
SOFT, the most reddened stars in or near the Ay maximum in the 
SOFT image are not detected in the SIRIUS image. 

The maximum visual extinction in the head computed for the 
SOFI data with the Mathis (1990) extinction law differs from the 
results obtained in Paper I. This is due to the maps using different 
binning, which causes the stars in the Nl frame to fall in different 
pixels. Because the Ay value in a single pixel is computed as an 
average of all the extinctions from the stars in the area of that 
pixel, the extinction values in the map can change when the pixel 
centre coordinates are changed. The Ay peak in the CG 1 head 
is sharp, which makes it especially sensitive to a change in the 
pixel coordinates. 



4.1.2. CG2 

The SIRIUS and SOFT Ay maps of CG 2 are shown in the upper 

and lower panels of Fig. 5, respectively. The lowest contour lev- 
els in the SIRIUS and SOFI maps are r."0 and 2'."0, respectively. 
The contour step is V^O. The observed distribution of extinction 
in CG 2 is similar to the one observed in CG 1. A sharp maxi- 
mum is seen in the globule head and a weaker maximum lies in 
the tail. The tail maximum is, however, relatively weaker when 
compared to the head clump than is the case in CG 1. In the 
SOFI map the visual extinction follows the filamentary structure 
that extends from the apex of the globule to the two bright stars 
in the NW direction. The Ay in the head peaks at 5^0 in SIRIUS 
and at 9^6 in SOFI. In the tail a maximum of 2"5 is observed. 
Outside the globule the Ay varies between 1™1 and 1"4. 




20=" 7'^16'^00= 40^ 20 




RA 02000) 

Fig. 5. Visual extinction map of CG 2. Upper panel: SIRIUS data. The 
contours are from 1™ to 5™ in steps of 1™. Lower panel: SOFI data. The 
contours are from 2" to 9™ in steps of 1". The grey scale is the same in 
both panels. 



4.2. Cloud mass and mass distribution 

The mass of a cloud can be estimated using the Bohhn et al. 
(1978) relation between the H2 column density and visual ex- 
tinction, N(H2) = 0.94 * lO^i cm-^ mag-\ where E(B-V)=3.1 is 
assumed. The extinction depends linearly on the molecular hy- 
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drogen column density, and thus Ay maps trace the cloud mass 
distribution. The choice of a reddening law and the adopted dis- 
tance to the object will influence the total mass of the cloud, but 
they do not influence the relative mass distribution within the 
Ay map. Thanks to its better resolution, SOFI data is used to es- 
timate the clump masses. However, the masses estimated from 
SIRIUS data do not diff'er significantly when using Ay contours 
that define an area equal in size. 

4.2.1. CG 1 

The CG 1 total mass as derived from the SOFI map, assuming 
a distance of 300 pc and using the BB88 extinction law, is 41.9 
Mq. Using the Rieke & Lebofsky (1985) reddening law the mass 
would be 45.7 Mq. The mass estimate for the Head clump inside 
the SOFI Av contour of 4.0™ is 6.9 Mq. In the Middle clump, a 
contour of 4.0*" yields 5.9 M©. The Tail clump mass inside the 
SIRIUS 2.5'" contour is 2.2 Mq. 

To investigate the MHead/^Toi mass ratio, the cloud is di- 
vided into two along the right ascension at 7^18™55^ Using 2.0™ 
as the lower Umit, the SOFI data gives the mass in the head as 
16.8 Mq and the total mass of the imaged portion of CG 1 as 
40.6 Mq. The same magnitude limit for SIRIUS yields 11.3 Mq 
for the head and 37. 1 Mg for the entire imaged area of CG 1 . The 
part of the tail that was imaged only by SIRIUS accounts for a 
mass of 8.6 Mq. 

SOFI data provides a mass fraction of 0.41 and SIRIUS data 
0.31 for the head. The discrepancy can be explained through the 
"extra" SIRIUS frame 14. Even though there is only a single, 
weak clump in this frame, there is a noticeable portion of mass 
spread uniformly throughout the frame. If the mass from the "ex- 
tra" SIRIUS frame is ignored, the mass fraction for SIRIUS >1'" 
becomes 0.34. The mass estimates may be uncertain by 50-100% 
owing to systematic errors in the Av values and the distance es- 
timate. However, the mass ratio between the head and the tail of 
the globule contains significantly less uncertainty because it is a 
differential measurement, and the error in the absolute value of 
the mass is likely to be fairly uniform throughout the cloud. 

Harju et al. ( 1 990) calculated a total molecular mass of 20- 
45Mo for CG 1 based on die '^CO -H2 and Ico-NCHi) conver- 
sion ratios. From the ^^CO map they determine that 25 % of the 
mass is located in the head. This agrees with our results for the 
mass fraction of the head, particularly when examining the larger 
SIRIUS-imaged area. 

4.2.2. CG 2 

The CG 2 total mass inside the 1™ contour is 31.0 M© for the 
larger SIRIUS Ay map and 13.1 Mq for the SOFI map. The SOFI 
contour 4™ gives a mass of 4.5 M© for the head. Dividing CG 2 
into two along the right ascension at 7*^1 5"45^ and using the 1™ 
contour, the mass fraction of the head is 0.62. 

The IRAS HIRES 12, 25, 60, and 100 //m images of CG 2 
are shown in the online Fig. B.6. The CG 2 total length in the 
IRAS HIRES images is 28', most of which is located outside our 
NIR frames. There is a local maximum in the tail in the IRAS 
25 fim image at 7''15"20'* but no corresponding feature in the 
SIRIUS Ay map (Fig. 5, upper panel). The IRAS 60 /im map has 
a local maximum ~5' from the head, which corresponds well to 
the local maximum in the SIRIUS Ay map. The portion of the 
tail not imaged in NIR contains one smaller clump, which has 
IRAS 60 and 100 fim fluxes less than those in the head. Including 



this would decrease the mass fraction from the value calculated 
above. 

4.3. Star formation in CG 1 and CG 2 

So far two stars, the pre-main- sequence star NX Pup (Reipurth 
1983) and 2MASS J07192185-4434551 which is embedded in 
the dense globule head 70" west of NX Pup (Paper 1), are known 
to have formed in CG 1. NX Pup is actually a close binary, and 
a third T Tau star lies just 7" NE of it (Bernacca et al. 1993; 
Brandner et al. 1995). The IRAS point source IRAS 07178-4429 
has earlier been often associated with NX Pup. However, its 
nominal position lies approximately halfway between NX Pup 
and 2MASS J07192185-4434551. As shown in Paper I, only the 
12 and 25 fim emission originates in NX Pup and the 60 and 100 
fim emission in the embedded 2MASS star. 2MASS J07192185- 
4434551 is referred to as CGI IRS 1 in the foflowing. 

Signs of further star formation in CG 1 and CG 2 were 
searched for by combining the NIR imaging and the photometry 
presented in this paper and images from the Wide-field Infrared 
Survey Explorer (WISE). WISE imaged the sky in four bands: 

3.4, 4.6, 12, and 22 fim (Wright et al. 2010). CG 1 and CG 2 
All-Sky Release Catalog data and WISE AU-Sky Release im- 
ages, public since March 2012, were retrieved via IRSA. 

The WISE images showed counterparts to three NIR excess 
objects, two in CG 1 and one in CG 2. None of these have opti- 
cal counterparts in either the DSS or the ESO and AAO R band 
images. These new objects are referred to as CG 1 IRS 2, CG 1 
IRS 3, and CG 2 IRS 1 in the foUowing. The SOFT and WISE 
photometry and archival catalogue data (Akari, IRAS, Spitzer 
if available) retrieved via IRSA for the five stars in CG 1 and 
CG 2 are listed in Table 1. The catalogued flux values were 
adopted when available, and for Spitzer the fluxes were extracted 
with MOPEX^ as aperture photometry. WISE Catalog data were 
expressed in magnitudes and converted into Janskys with the 
zero-point magnitudes listed in the WISE All-Sky Data Release 
Product Explanatory Supplement^. 

4.3.1. NX Puppis 

Schoeller et al. (1996) derived the masses » 2.0 Mq and 1.6-1.9 
Mq and ages 3 - 5 * 10^ and 2-6*10^ years for the binary NX 
Pup components A and B, respectively. Brandner et al. (1995) 
give component B a wider range of mass, 1.2 - 2.5 M© and a 
possibility to be as young as 0.3 * 10^ yr. The third component 
7'.'0 away from the binary, NX Pup C, has a mass of 0.30 Mq. 
We have not modelled NX Pup as it has already been done by 
Schoeller et al. (1996) and Brandner et al. (1995) and the binary 
nature of the system is unresolved in our observations. However 
we note that the IRAS 60 and 100 jum fluxes deviate strongly 
from the modelled binary spectral energy distribution (SED) in 
Brandner et al. (1995). This is expected because the mid and far 
IR signal is not cormected with NX Pup at all. 

4.3.2. CG 1 IRS 1 (2MASSJ071 921 85-4434551) 

This object was discussed in Paper 1 where it was identified as a 
young stellar object. A molecular hydrogen object MHO 1411, 
which is a probable obscured Herbig-Haro (HH) object, Ues ~ 
90" west of CG 1 IRS 1. JHKs photomefiy of CGI IRS 1 is not 

' http://irsa.ipac.caltech.edu/data/SPITZER/docs/dataanalysistools/ 
tools/mopex/ 

^ http://wise2.ipac.caltech.edu/docs/release/allsky/expsup/ 
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possible as it is totally obscured in J and H bands and only re- 
flected light is seen. The central object is visible in Ks but even 
here most of the light is due to scattering in the circumstellar ma- 
terial. Hence we have excluded these bands from our photome- 
try. The IRAS Point Source Catalog fluxes of IRAS 07178-4429 
at 60 and 100 yum are attributed to CGI IRS 1 in Table 1. Akari 
observations are both from the FIS Bright Source Catalogue and 
IRC Point Source Catalogue. 



4.3.3. CG 1 IRS 2 (WISE J071 846.52-443252.9) 

This obj ect is located in the tail part of CG I , in what is identified 
as an area of extended surface brightness in Fig. 1, and at the NE 
edge of the SOFI Middle clump in Fig. 4. The colours derived 
from SOFI observations are 7 - // = 1 .74 and H - Ks, - \ .24 
which put it in the upper right-hand corner of the J-H, H-Ks 
colour-colour diagram in the left panel of Fig. 3. The WISE im- 
ages show a matching object in all wavelengths. 



4.3.4. CG 1 IRS 3 (WISE J071 81 4.71 -443420.0) 

This object is detected in the tail of CG 1. It is not included in 
the final SOFI photometry catalogue because SE classified it as 
a non-stellar object. WISE also flags it as an extended object, 
possibly due to confusion noise, since there are several stars lo- 
cated within ~7" of IRS 3 in the NIR images. The highest res- 
olution of WISE is 6'.'1 at 3.4 /im. In fliis band WISE shows a 
bright peak at IRS 3 and at one of these close stars, but the WISE 
Source Catalog lists only IRS 3. We have included CG 1 IRS 3 
in our NIR excess object list because of its location in CG 1 and 
because it can be seen in the WISE 22 fim image. 



4.3.5. CG 2 IRS 1 (WISE J071 604.50-435742.6) 

This star lies in the head of CG 2 between the globule apex 
and the region of maximum visual extinction. SOFI photome- 
try gives the NIR colours J - H - 1.73, H - Ks - 1.45. In the 
JsHKs images it is clearly red and appears slightly extended. In 
the 22 f^m WISE image, this is the brightest object in the CG 2 
head. WISE 12 /xm image of CG 2 is shown in Fig. 6. Contours 
of the HIRES-processed 100 fim map and of visual extinction 
derived from the SOFI data are also shown. The large artefact 
in the CG 2 tail is due to a reflection from a bright nearby star. 
The IRAS 25 i^m grey scale and IRAS HIRES-processed 60 //m 
and 100 /urn contours are shown in the online Fig. B.7. Bhatt 
(1993) associated IRAS 07144-4352 with the head of CG 2. 
Even though only low-quahty IRAS PSC 25 jum flux is available 
for the IRAS source, a point source is visible in the HIRES 25 
fim grey scale image. The position of this point source coincides 
wifli that of CG 2 IRS 1. The HIRES 100 //m contours are cen- 
tred SW of the 25 fim source and trace the dense globule core 
(Fig. 6) rather than IRS 1. Similar to IRAS 07178-4429 in flie 
CG 1 head, which traces emission both from NX Pup and CG 1 
IRS 1 (see Paper I), IRAS 07144-4352 also consists of two com- 
ponents: the molecular core and CG 2 IRS 1. The IRAS 60 fim 
flux traces, at least partly, CG 2 IRS 1 and the 100 yum flux the 
dense globule core. The Akari FIS catalogue contains a source 
that has coordinates that are offset ~29" southeast of the WISE 
coordinates. Neither the IRAS nor the AKARI far infrared fluxes 
are associated with CG 2 IRS 1 in Table 1. 



4.3.6. SED modelling 

A SED fitting tool by Robitaille et al. (2007) is available online. 
The SOFI, WISE, IRAS, and Spitzer data in Table 1 were used 
as the input. The errors were set to 10% although the true errors 
are smaller than that. The results are only indicative especially 
because of the lack of FIR and submm data. Only CG 1 IRS 1 
has a measured flux at 100 fim, but for the rest of the sources the 
longest observed wavelength is the WISE 22 fim. The distance to 
the objects was allowed to vary between 300 and 450 pc to reflect 
the literature estimates for CG 1. For CG 1 IRS 3, a range of 
200-600 pc was allowed because of its location far in the CG 1 
tail. The interstellar extinction was allowed to vary from 0.1™ 
to a maximum estimated by using the Ay maps. The age, mass, 
temperature, and luminosity of the central source. Ay, distance 
and the;^^^ for the best fit are listed in Table 2. The mass and age 
ranges for all good fits (i.e. fits that have^^ ~xlest ^ ^ P^'^ ^^'^^ 
data point) are also listed. The SED plots with the fits are shown 
in Fig. 7. The uncertainty of the long wavelength part of the SED 
due to the lack of data points is visible in the plots. Online Fig. 
B.8 shows the age-mass relation derived from the good fits. All 
fits point to low-mass stars. 

The SED for IRS 1 in the head of CG 1 has a shape of a class 
I YSO. The best fit gives a mass of 2.8 and an age of 0.5 
Myr, which would make IRS 1 more massive than the individual 
NX Pup components. The other good fits have masses below 1 
Mq, and their ages are in the range of 3 - 217 kyr. The age range 
covers two orders of magnitude, but even with the highest age 
estimate IRS 1 is a second-generation star compared to NX Pup. 
The star formation efficiency in the head of CG 1 is about 25% 
when using the Schoeller et al. (1996) mass estimates for the AB 
components of NX Pup, the SOFI Mnead estimate and the mean 
IRS 1 mass from flie SED fits (0.55 M©). 

The NIR excess objects CG 1 IRS 2 and IRS 3 have similar 
SEDs and are best fitted with YSOs with ages in the range ~l-7 
Myr and masses in the range ~ 0.2-1.3 M©. The age is similar 
to the age of NX Puppis. 

In CG 2, the stellar age in the fits ranges from 1.2 kyr to 0.14 
Myr and the masses from 0. 1 1 to 0.29 Mq. The age of this object 
is comparable to the age of IRS 1 . The lack of a large amount of 
obscuring circumstellar material in CG 2 IRS 1 suggests that this 
object is more evolved than CG 1 IRS 1. 

4.4. Morphology 

In addition to the new J?HKs observations, satellite data (IRAS, 
Spitzer, WISE) and molecular line data from Paper 1 and Harju 
et al. (1990) are used to study the structure of CG 1. Only NIR, 
IRAS, and WISE data are available for CG 2. The large-scale 
NIR structure of CG 1 and CG 2 has been discussed in Sects. 
3.1.1 and 4.1.1. 

Large polycycUc aromatic hydrocarbon (PAH) molecules 
emit strong emission lines at 3.3, 6.2, 7.7, 8.6, 11.2, 12.7, and 
16.4 fim (Tielens 2008). PAH emission is thus observed in 
Spitzer IRAC bands 3.6, 5.8, and 8.0 fim, and in WISE bands 
3.4 and 12 fim. The amount of PAH emission in the diflxise in- 
terstellar medium observed in the Spitzer IRAC bands has been 
investigated by Flagey et al. (2006). Very small grains (VSGs) 
emit thermally in the mid infrared (MIR) and are expected to be 
seen only in the Spitzer MIPS 24 and 70 fim and IRAS 25, 60, 
and 100 fim filters. Thermal emission from the large, cold grains 
is expected only in the IRAS 100 fim filter but even here a large 
fraction of the observed emission is due to small grains. Also 
molecular hydrogen has numerous rotational fines in the near- 
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Fig. 6. WISE 12/im image of CG 2. Arbitrary contours of IRAS 100 jim HIRES processed image (black lines) and the extinction derived from the 
SOFT data (grey lines) are also shown. The strong artefact in the CG 2 tail is a reflection of a nearby bright star. 

Table 1. The photometry for NX Pup and the new NIR excess candidates in CG 1 and CG 2. 





NX Pup 


CG 1 IRS 1 


CG 1 IRS 2 


CG 1 IRS 3 


CG 2 IRS 1 




590 (8.579) 




0.0235 (19.58) 


0.0435 (18.91) 


0.25 (17.01) 


w 


1250 (7.825) 




0.0749 (17.84) 


0.0943 (17.59) 


0.791 (15.28) 


K" 


2470 (6.080) 




0.153 (16.60) 


0.158 (16.56) 


1.96 (13.83) 


3.4* 


5184.62 (4.440) 


7.75 (11.504) 


0.281 (15.104) 


0.453 (14.586) 


6.923(11.626) 


3.6' 




13.94 








4.5' 




23.39 








4.6* 


8593.90 (3.252) 


23.05 (9.681) 


0.369 (14.170) 


0.579 (13.680) 


12.853 (10.315) 


5.8^ 




22.29 








8.0" 




34.24 








9.0" 


5197 










12* 


5047.81 (1.994) 


144.91 (5.849) 


2.403 (10.300) 


2.271 (10.361) 


24.769 (7.767) 


12-* 


6682 








250* 


18^ 


4652 


1029 








22* 


4816.84 (0.599) 


2236 (1.432) 


8.279 (7.511) 


11.261 (7.177) 


58.879 (5.381) 


24-^ 




2092 








25'' 


7604 








370'' 


60'' 




13120 






440'' 


65« 




7697 












9314 








9W 




9672 






656 


100^ 




33590 






8880'' 


140? 




17370 






6259 


160* 




15570 









Notes. The flux unit is mjy. For JHK and WISE we also Ust flie magnitude in parenthesis. SOFI. <*> WISE. Spitzer IRAC. IRAS. <"* Akari 
IRC. <^ Spitzer MIPS. <«' Akari FIS. Data for the offset IRAS source 07144-4352, flagged U (upper limit) for 12 and 25 jum, B and A for 60 and 
100//m, respectively. Used to check the SED model behaviour. 



Table 2. The results from the best stellar model fits, where mass and age ranges include all models with^ ~Xi,est ^ P^"" '^^'^ point. 





CG 1 IRS 1 


CG 1 IRS 2 


CG 1 IRS 3 


CG 2 IRS 1 


Mass (Mo) 


2.84 


0.50 


0.21 


0.14 


Mass (Mo) 


0.24-2.84 


0.22-1.22 


0.10-1.49 


0.11-0.29 


Age (yr) 


4.85(5) 


2.15(6) 


4.11(6) 


1.37(5) 


Age (yr) 


2.8kyr-485kyr 


1.9-7.2Myr 


0.1-9.5Myr 


1.2kyr-140kyr 


Temperature (K) 


4774 


3759 


3234 


2962 


hroT (Lo) 


1.62(1) 


4.05(-l) 


8.98(-2) 


3.65(-l) 


Av 


10.00 


3.59 


4.25 


6.70 


Distance (pc) 


331 


331 


209 


302 


Xtot 


41.1 


33.6 


10.7 


9.6 
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Fig. 7. SED fits for our NIR excess sources. NX Pup is omitted as it lias been previously modelled (see text). The filled circles are the observed 
fluxes and the triangles mark the fluxes where there is only an upper limit. The solid line marks the best fit, and the grey lines shows fits with 
- xle^, < 3 per each data point. The dashed line indicates the unreddened SED of the stellar photosphere of central source in the best fitting 
model. 



and mid-infrared. The high spatial resolution of the Spitzer and 
WISE images allows us to study the detailed structure of CG 1 
and CG 2 and especially to estimate whether the structures are 
due to PAH emission or not. 



cones, which are best seen in the IRAC 3.4, 4.5, and 8.0 /urn 
colour composite image shown in online Fig. B.13. The cones 
open to east and SW from CG 1 IRS 1. The opening angle of the 
cone to the east is larger than that of the one opening to SW. 



4.4.1. CG 1 

The SOFI JsHKs and H2 (1-0) S(l) emission structure of the 
CG 1 head was discussed in Paper I. The individual Spitzer 
IRAC images at 3.6, 4.5, 5.8, and 8.0 /mi and MIPS images at 
24 and 70 //m are shown in the online Figs. B.9 and B.IO, re- 
spectively. The Spitzer images cover only the head of CG 1 . The 
WISE images of CG 1 at 3.4, 4.6, 12, and 22 jjm are shown in 
the onhne Fig. B.ll. A false-colour WISE image is shown in 
Fig. 8. The location of the most notable features seen in the head 
of CG 1 in WISE and Spitzer are marked in the online Fig. B.12. 
Two cone-like structures symmetrically around CG 1 IRS 1 are 
seen in the IRAC images, with a bright streamer that originates in 
CG 1 IRS 1 and falls on the axis of the SW-pointing cone. This 
feature is seen best at wavelengths longer than 5.8 fim. Below 
NX Puppis is a bright reflection nebula that is seen only at 3.6 
and 70 /im. Surface brightness maxima are observed in the tip 
of the SW streamer and NW of CG 1 IRS 1 and are seen in the 
MIPS 24 yum image. The NW patch is marked with an elUpse. A 
soUd white circle marks MHO 1411. 

The cones: The spatial resolution of the WISE images is 
insufficient to resolve the thin filaments originating from CG 1 
IRS 1 seen in the SOFI JsHKs images (Paper I, Fig. 3). However, 
they are faintly visible in the IRAC 3.4 and 4.5 /im images where 
another SW pointing filament is seen. These filaments form two 



WISEJSW streamer: The bright, large-scale streamer 
reaching SW from CG 1 IRS 1 inside the SW cone will be called 
W1SE_SW. This feature is strong in the WISE 12 and 22 and 
IRAC 8.0 and MIPS 24 ij.m images. It is also faintly visible at 
3.4 and clearly detected in the 5.8 fim image. The maximum 
WISE_SW surface emission seen in the 24 /im image in the tip 
of the streamer lies just outside the other Spitzer images but is 
detected in the WISE 3.4, 12, and 22 /mi images. It is also seen 
in the HIRES-enhanced IRAS 60 fim image shown in the onUne 
Fig. B.14 (bottom left). Contour plots of WISE 3.4 and 22 //m 
images overlaid on the SOFI 7s image are shown in the online 
Fig. B.15. The contours dehneating WISE_SW in Fig. B.15 are 
above and clearly offset from the filament seen in 7s. WISE_SW 
and the bright filament seen in the JHKs imaging are thus dif- 
ferent features. The surface brightness maximum in the west- 
ern tip of WISE_SW has a corresponding patch on the edge of 
ffie globule head NW from CG 1 IRS 1. The NW patch has a 
lower surface brightness but is visible in the same images as the 
WISE_SW streamer except for MIPS 24 /im where the patch is 
outside the image. It is suggested that the WISE_SW streamer 
and the patch are due to PAH emission in the WISE 12 fira band. 
The emission seen in the IRAS HIRES 60 /im image would be 
due to thermal emission from very small grains. The source of 
excitation of the PAH emission and the heating mechanism of 
the VSGs is not known. 
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MHO 1411: This faint molecular hydrogen object was first 
reported in Paper I. A faint object is also seen in the IRAC 3.6 
and 4.5 jum and in the WISE 4.6 jum image. Its existence has thus 
been independently confirmed. 

Reflection nebulosity below NX Pup: The bright surface 
emission coincides with a small, localized component in the 
CG 1 C'^O mapping reported in Paper I where the feature was 
referred to as C^^O_SE. The nebulosity can be seen in the IRAC 
3.6 jim image (online Fig. B.9, top left) and the MIPS 70 fim im- 
age (online Fig. B.IO, right panel). The IRAC data strengthens 
the claim the nebulosity in 3.6 fim is due to NIR scattering of ra- 
diation coming from NX Pup. The surface brightness observed 
at 70 yum is caused by the thermal emission of very small grains. 

CG 1 Tail: The large-scale structure of the tail as seen in the 
NIR extinction and CO emission (Harju et al. 1990) agrees well 
with the tail seen in the optical images. A new structure is seen 
in the WISE 3.4, 12, and 22 yum images (online Fig. B.l 1) and in 
the HIRES -enhanced IRAS maps (online Fig. B.14). Extended 
surface brightness at the southern edge of the tail starts about 
13' west of the head and extends ~10' west. This WISE surface 
brightness is not seen in the NIR images. In the following this 
feature is referred to as TaiLSouth. The infrared excess star CG 1 
IRS 3 lies at the very eastern edge of TaiLSouth. 

TaiLSouth resembles a cometary globule superposed on 
CG 1. However, the lack of both observed extinction and of 
strong CO emission (Harju et al. 1990) suggests that it is not 
a second cometary globule. More likely TaiLSouth is emission 
from PAHs and/or VSGs. IRAS HIRES data shows a local max- 
imum at the position of TaiLSouth at 12, 25, and 60 /urn. In the 
WISE images TaiLSouth is strongest at 12 /xm and weakest at 
3.4 jum. These two wavelengths could be explained with PAH 
emission, but the emission observed at 22 fim and in the IRAS 
HIRES 25 and 60 jum images requires VSGs. This raises an- 
other question: What is the origin of excitation of the PAHs and 
the heating for the VSGs? It cannot be the general interstellar 
radiation field (ISRF) from the Galactic plane as TaiLSouth lies 
on the shadow side of the CG 1 tail with respect to the Galactic 
plane. In addition, the head of CG 1 shadows TaiLSouth from the 
emission from NX Pup and the O star ^ Pup in the centre of the 
Gum Nebula. However, the multiple system Vel in the central 
part of the Gum Nebula is located so that a line drawn from the 
system through the tip of WISE_S W would follow the NE edge 
of TaiLSouth. In this view the SW streamer shields the part of 
the CG 1 tail that is closest to the head and leaves the rest of the 
tail exposed to the UV radiation of Vel. In the tail, the IRAS 
12 and 25 fim emission peak is separated from the emission in 60 
and 100 fim. Similar behaviour has been detected e.g. in L1780 
by Ridderstad et al. (2006). They suggest that this indicates two 
spatially segregated dust populations. 

4.4.2. Does CG 1 IRS 1 drive an outflow? 

CG 1 IRS 1 has the appearance of a Class I YSO (see Fig. 7 and 
Paper 1). Such objects typically drive molecular outflows and 
jets. However, in this case there is very little evidence of an out- 
flow in the available molecular line data presented in Harju et al. 
(1990). An indication of a very modest outflow wing is seen in 
the '^CO (3-2) spectra covering the immediate surroundings of 
CG 1 IRS 1 (Paper I). The only firm indicator of an outflow is 
MHO 1411, which lies about 90" west of CG 1 IRS 1 . In addition 
to the molecular hydrogen object, the JHKs and Spitzer IRAC 
3.6 and 4.5 fim images reveal two cones that could delineate the 
walls of a cavity carved by an outflow. The cone edges are vis- 
ible in the IRAC 4.5 fim image and are thus not due to PAH 



emission, leaving scattering as the viable option. IRAC observa- 
tions around HH 46/47 presented in Velusamy et al. (2007) show 
that the outflow cavity edges can be visible in the IRAC images. 
The WISE_SW streamer is puzzling because it lies between the 
possible cavity walls. At 12 /im it could be due to strong PAH 
emission but the detection at 22 and 60 jum suggests that it is 
also associated with warm VSGs. Considering the new informa- 
tion provided by Spitzer imaging, it is suggested that CG 1 IRS 1 
is driving an outflow, and the outflow cavity walls are seen in 
reflected light. However, the nature of WISE.SW and the very 
weak CO outflow need to be investigated further in detail. 

4.4.3. CG 2 

CG 2 WISE images at 3.4, 4.6, 12, and 22 /im are shown in on- 
line Fig. B.16. The artefact in the centre of all images is caused 
by a very bright field star south of the field. The globule head 
is faintly detected at 3.4 ^um and hardly visible at 4.6 jum. At 12 
and 22 fim the surface brightness is stronger and coincides ap- 
proximately with the extended surface brightness of the SIRIUS 
and SOFI images (Figs. B.3 and 2). The NIR excess object CG 2 
IRS 1 is the brightest object in the CG 2 head at 22 jum. The nar- 
row tail is mostly covered by the artefact, but a lane of emission 
is also seen both to east and west of it. This emission lane can be 
seen in IRAS HIRES contours (Fig. B.6) at 25, 60, and 100 //m. 
The WISE artefact covers the western maximum emission at 60 
and 100 fxm (Fig. 6). The total length of CG 2, from the head to 
the end of the emission lane is ~28'. 

No information on the extinction in the tail is available be- 
cause it is not included in the NIR observations. The tail is sepa- 
rated into individual clumps in the HIRES images, but this may 
possibly be because the HIRES maximum entropy algorithm 
tends to form compact structures. At the moment it is not possi- 
ble to say whether the MIR and FIR emission seen in the globule 
tail is only due to warm small dust particles or if it is because 
of a combination of both very small and larger sized particles. 
The extension of the NIR imaging, molecular line observations 
or submm dust continuum observations are needed to solve this 
question. 

4.5. CG formation mechanism 

The mass distribution in the shocked cloud is discussed in the 
literature mostly on a qualitative level, but some estimates have 
been given. Heathcote & Brand (1983) find that in the SN shock 
case the fraction F of mass that goes into the tail depends on the 
velocity of the flow of the gas behind the external shock. For a 
typical case they find F ~ 0.5. In the RDI simulations of Lefloch 
& Lazarefi" (1994), 80% of the mass remains in the head after 
the initial collapse phase. The cometary phase begins after the 
collapse phase and lasts 90% of the estimated 3 Myr lifetime 
of the CG. The flux of the ionizing radiation can also affect the 
mass distribution in RDI (Gritschneder et al. 2009; Bisbas et al. 
2011). 

The mass fraction of the CG 1 head has a lower limit of ~0.4 
from the SOFI observations and ~0.3 from the SIRIUS. Star for- 
mation in CG 1 has taken place first in the globule head with 
NX Puppis as the first generation and CG 1 IRS 1 the second, 
triggered generation. The position of these stars near the glob- 
ule's front edge is in accordance with the RDI-induced forma- 
tion (Sugitani et al. 1989; Lee & Chen 2007; Bisbas et al. 201 1). 
Taking also the masses of these stars into account, the fraction 
of mass in the CG 1 head is raised by about 7-10%. As the star 
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Fig. 8. WISE CG 1 false-colour image. The 4.6, 12, and 22 jim bands are coded in blue, green, and red, respectively. 



formation efficiency (SFE) is typically of the order ~few % in 
dark clouds, the original mass of CG 1 head must have been 
higher before star formation took place. This pushes the mass 
fraction towards the values from the RDI model. Even though 
the CG 1 head-to-tail mass ratio may agree with the RDI for- 
mation scenario, the disrupting effects of star formation prevent 
firm conclusions. The globule is in too evolved a stage to reflect 
the original mass distribution after the triggering event. 

The observed head-to-total mass ratio in CG 2 is 0.6. This 
ratio must have been higher before the formation of CG 2 IRS 1 
in the head. This would be a strong argument in favour of RDI 
being the globule formation mechanism. However, the observed 
area does not cover the whole tail. The appearance of the tail 
is by far not as extended as in CG 1. Further observations of 
the tail are needed before a definitive value can be given for the 
head-to-tail mass ratio in CG 2. 



5. Summary and conclusions 

We have analysed NIR imaging and photometry and available 
archival data to study the detailed structure and star formation 
in the cometary globules CG 1 and CG 2 located in the Gum 
Nebula. The NIR observations were done in the NIR /, H, and 
Ks bands with SOFI at the NTT telescope and SIRIUS at IRSF 
In addition, data from IRAS, WISE, Spitzer, and AKARI tele- 
scopes were retrieved from online archives. 

- We find two new NIR excess objects in CG 1, and one in 
CG 2. CG 1 IRS 2 and CG 1 IRS 3 are located in the tail, 
whereas the previously known young stellar CG 1 IRS 1 is 
located in the head of the globule (Paper I). The formation of 
the pre-main sequence star NX Pup in CG 1 is already known 
(Reipurth 1983). The object CG 2 IRS 1 located in the head 
of the globule is the first indication of star formation taking 
place in CG 2. 

- According to SED fitting, CG 1 IRS 1 is a class I source. The 
SED fits to the rest of the NIR excess sources indicate that 
they are also young, low-mass stars. The spread of the fitted 
ages for each NIR excess object is, however, large. This is 



largely due to lack of FIR and submm photometry data for 
all objects except CG 1 IRS 1. 

- CG 1 IRS 1 in the head of CG 1 is a second-generation star 
and is indicative of triggered star formation. 

- CG 1 IRS 1 is likely to drive a bipolar outflow. Structures 
similar to outflow cones are detected in the imaging. The 
wider eastern lobe opens towards NX Puppis and the western 
lobe opens to SW. Mysteriously, no strong outflow features 
are seen in the available molecular data, but a bright streamer 
is seen on the axis of the SW cone in MIR. 

- WISE images and IRAS HIRES-processed images suggest 
the presence of excited PAH particles and warm, very small 
grains in the CG 1 tail. 

- CG 1 is in far too evolved a stage to reflect the globule mass 
distribution at the time of its formation. Thus no firm conclu- 
sion on the original tail-to-head mass ratio can be obtained. 
However, when considering the star formation history in the 
globule, the mass distribution of CG 1 is indicative of RDI 
triggered formation. 

- The observed tail-to-head mass ratio in CG 2 is in accor- 
dance with the RDI triggered formation. However, since the 
observations do not cover the tail in whole, further observa- 
tions are needed before a final conclusion can be drawn. 
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Appendix A: Data reduction 

The SOFI and SIRIUS data were reduced in the same man- 
ner as in Paper I. IRAF^ and the external XDIMSUM package 
was utilised. For SOFT, each frame was corrected for cross talk. 
Then bad pixels were masked and cosmic rays removed in all 
frames. Sky-subtraction was done using the two temporally clos- 
est frames. An object mask was created for each frame to mask 
the field stars. Running sky-subtraction with the object masks 
also created hole masks for each frame. The flat field and illumi- 
nation correction files available at the ESO SOFI website were 
used on the SOFI sky-subtracted images. Evening sky flats were 
used to flat-field the SIRIUS data. The bad pixel mask, cosmic 
ray masks, and hole masks were combined to create rejection 
masks, which were used when co-adding the images. The coor- 
dinates were tied to the 2MASS Point Source Catalog (PSC). 

For each colour the objects were extracted from the reg- 
istered frames using Source Extractor v. 2.5.0 (SE, Bertin & 
Amouts 1996). SE fits an aperture to each object, and this is 
used to extract the flux information. The SIRIUS photometry 
was calibrated using 2MASS PSC. The SE photometry of each 
colour and field was caUbrated using 10-14 2MASS PSC stars 
per field by linearly scaling the frames to the PSC star level. 
The difference in the magnitudes of the observed field and the 
PSC was < ±0.01™ in all colours. SOFI-magnitude zero points 
were fixed using standard star observations. The SOFI magni- 
tudes extracted by SE were converted into the Persson system 
(1998) and from there to 2MASS magnitudes, as described by 
Ascenso et al. (2007). The SOFI magnitude scale was checked 
using stars in common with the 2MASS PSC. 

Several experimentally found selection rules were used to 
discard galaxies and stars with poor photometric data from the 
catalogues produced by SE. SE photometry of SIRIUS images 
were filtered as follows. Objects were removed if they were 
closer than 80 pixels to the frame border, SE flagged their flux 
data, their elongation was more than 1.3 or if any colour had a 
photometric error calculated by SE of larger than 0^15. SE also 
provides keyword CLASS, which assigns a number between 0.0 
and 1.0 to describe how star-like the object is (1.0 for stars, 0.0 
for galaxies) in each colour. This was utilised so that objects 
that had a summed-up star index less than 1.6 were discarded. 
The selection rules to exclude non- stellar objects and stars with 
doubtful photometry in the SOFI SE photometry were the ob- 
jects closer than 30 pixels (0 for CG 2) of the edge of the reg- 
istered frame, those SE flagged to have bad flux data, objects 
with an elongation larger than 1.3 (1.25 for CG 2), or a magni- 
tude error over 07'15 in any of the three bands. Objects whose 
summed-up star indices were less than 2.70 were left out. After 
a visual inspection, any objects with two unresolved stars in the 
extraction aperture were removed. The OFF field SE photome- 
try was filtered the same way. In practice the formal error of the 
remaining objects is less than 0^1 in aU colours. 



Appendix B: Online data 



' IRAF is distributed by the National Optical Astronomy 
Observatories, which arc operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 
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Fig. B.l. Colour-coded SIRIUS image of CG 1. The J, H, and A"s bands are coded in blue, green, and red. 




Fig. B.2. SOFI Js, H, and Ks band images of CG 1 . 
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Fig. B.3. Colour-coded SIRIUS images of CG 2. Tlie J, H, and A^s bands are coded in blue, green, and red. The object CG 2 IRS 1 is indicated 
with a box. 
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Fig. B.4. Sirius J, H, and Ks images of CG 2. 
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Fig. B.5. SOFT 7s, H, and ATs band images of CG 2. 
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Fig.B.6. IRAS HIRES processed images of CG 2 in 12, 25, 60, and 100 yum. 
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Fig. B.7. IRAS 25 jum HIRES processed image of CG 2. Overplotted are the IRAS HIRES processed 60 nm (black) and 100 nm (red) contours 
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Fig. B.8. Age and mass in the SED fits. The filled grey circles show the distribution of models in the model grid, and the black filled circles shows 
the distribution of the best-fitting models i.e. with^^ ~ xlest ^ ^ P^'^ ^^^^ '^^^^ point. 
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Fig. B.9. Spitzer IRAC 3.6, 4.5, 5.8, and 8.0 //m images of CG 1 . 
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Fig. B.IO. Spitzer MIPS 24 (left) and 70 jum (right) image of CG 1. 
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Fig. B.ll. WISE 3.4, 4.6, 12, and 22 //m images of CG 1 . 
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Fig. B.12. Observed structures in the CG 1 head over-plotted on the false-colour SIRIUS image. The thin, solid arcs starting from the YSO denote 
the possible outflow cones. The dashed semi-circle below the brightest objects NX Puppis marks the reflection nebula. The solid circle marks 
the location of the molecular hydrogen object MHO 1411. The ellipse in the top right-hand corner denotes the location of the patch of surface 
brightness, and the extended ridge (solid line) in the bottom right-hand corner marks the bright WISE_SW streamer filament seen in the Spitzer 
images. 
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Fig. B. 13. IRAC falsecolour image of CG 1. The 3.6 is coded in blue, 4.5 yum in green, and 8.0 yum in red. 
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Fig.B.14. IRAS fflRES processed images of CG 1. Top left: 12 pm. Top right: 25 fim. Bottom left: 60 yum. Bottom right: 100 yum. 
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Fig. B. 15. 7s image of CG 1. Left: The WISE 3.4 yUm contours have been overplotted. The dashed lines indicate the strongest NX Pup diffraction 
spilces. Right: The WISE 22 yum contours have been overplotted. The contour levels are arbitrary. 




Fig. B. 16. WISE CG 2 in 3.4, 4.6, 12, and 22 //m images. The strong artefact in the tail is caused by a nearby bright star. 



